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Abstract 

We consider clean cylindrical nanostructures with an applied longitudinal 
static magnetic field. Without Coulomb interaction, the field induces, for 
particular values, points of degeneracy where a change of ground state takes 
place due to Aharonov-Bohm effect. The Coulomb potential introduces inter- 
action between the electronic configurations. As a consequence, when there is 
degeneracy, the ground state of the system becomes a many body state - un- 
able to be described by a mean-field theory - and a gap is opened. To study 
this problem, we propose a variational multireference wave function which 
goes beyond the Hartree-Fock approximation. Using this ansatz, in addition 
to the avoided crossing formation, two other effects of the electron-electron 
interaction are pointed out: (i) the long-range part of the Coulomb potential 
tends to shift the position in magnetic field of the (avoided) crossing points 
and, (ii) at the points of (near) degeneracy, the interaction can drive the sys- 
tem from a singlet to a triplet state inducing new real crossing points in the 
ground state energy curve as function of the field. Such crossings should ap- 
pear in various experiments as sudden changes in the response of the system 
(magnetoconductance, magnetopolarisability,...) when the magnetic field is 
tuned. 



I. INTRODUCTION 



It is known, since the early days of quantum mechanics, that the low energy electronic 
properties of aromatic molecules are very sensitive to a magnetic field applied perpendicu- 
larly to their planes [1]. The field breaks the time reversal symmetry and therefore induces 
an electronic current running around the circumference of the ring. This is the persistent 
current, arising due to Aharonov-Bohm effect [2]. This is an equilibrium phenomenon, 
periodic in magnetic flux with period 0o = hc/e, the flux quantum [3]. But, to be able 
to measure it requires systems with characteristic lengths in the nanoscopic or mesoscopic 
scale. Indeed, on the one hand, the motion of the electrons has to stay coherent which is 
possible for systems smaller than the electronic coherence length. But, on the other hand, 
to cover a full period of magnetic flux requires rings with sufficiently large diameters. This 
is definitely not the case for the usual aromatic molecules where field as large as 10 5 tesla 
are needed to observe the periodicity. However, nowadays, many man-made systems in the 
appropriate length scale, are available in various forms: isolated or ensemble of metallic 
or semiconducting rings [4] and, more recently, carbon nanotubes which are large organic 
cylinders [5,6] or, last, rings of carbon nanotubes [7,8]. Therefore, the study of persistent 
current have regained lot of interest during the last ten years. 

Rings, made of usual metal or semiconductor, are studied intensively since the nineties 
and a few experiments, motivated by early theoretical prediction [9] , have measured a sizable 
persistent current in some systems [4]. But, neither the magnitude of the current - one or two 
orders of magnitude larger than expected - nor the diamagnetic sign of the response measured 
experimentally can be explained by existing theories. Most of the theoretical efforts were 
devoted to the study of the interplay between disorder and Coulomb interaction, but without 
convincing conclusions up to now [10]. 

Carbon nanotubes were discovered by Iijima in 1991 [5]. They are fascinating materials 
whose electronic properties are determined in a unique way by the topology of their lattice: 
they are rolled up strip of graphite sheet whose, depending on their diameter and chirality, 
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can be either metallic or semiconducting [11]. In any case, the electronic spectrum of these 
systems seems to be very sensitive to an applied magnetic field suggesting large orbital mag- 
netic response [12,13]. Indeed, recently strong field effects have been seen in measurements 
of the tunnelling conductance of multiwall carbon nanotubes [6]. Magnetoconductance mea- 
surements of rings of carbon nanotubes have also been done [7]. 

In order to better understand these experimental results, additional studies of electronic 
multichannel systems in presence of Coulomb interaction are needed. This is clear for 
ballistic systems but, even for metallic rings which are in the diffusive regime, a proper 
understanding of the clean case should help to find a formalism able to treat disorder and 
interaction on equal footing. In this work, we consider cylindrical systems made of rolled 
square lattices (multichannel systems), without disorder but with Coulomb potential, long 
or short range. We are interested in how the Coulomb potential can affect the ground state 
properties of these systems when a static magnetic field parallel to the cylindrical axis is 
tuned, and hence, on the behaviours of the persistent current. Similar studies were done in 
the past but, either for pure ID systems [14] i.e. systems with only one electronic channel 
(rings), or for multichannel systems (cylinders) but with strong disorder using first order 
perturbation theory or the Hartree-Fock approximation [15]. The conclusions obtained here 
are not contained in these works. 

For cylinders, in the pure case and without electron-electron interaction, an axial mag- 
netic field induces many level crossings at different values of the magnetic flux [16]. At zero 
temperature, for a certain density of particles, the ground state energy shows also crossing 
points for particular values of the field where the ground state is degenerate. With the 
Coulomb interaction, there exist interaction between these degenerate states resulting in 
avoided crossings; in other words, at least at the vicinity of the crossing points the ground 
state becomes a true many-body state unable to be described by any mean field treatment. 
In this work, we propose a 'minimal' variational wave function to deal with this particular 
problem, going beyond a simple Hartree-Fock calculation. As a result, we find, in addition 
to the avoided crossing formation, two others effects caused by repulsive interaction, (i) The 
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positions of the (avoided) crossings points are shifted in magnetic field; this effect is due to 
the non-equal Hartree contributions of the different components of the interacting ground 
state, (ii) Because of the presence of degeneracy, the Hund's rule may drive the system to a 
triplet state: it follows sequences of singlet — > triplet — > singlet ground states which should 
be seen in various type of experiments such as measurements of the magnetoconductance, 
for instance. 

We believe the conclusions of this work rather general providing that the system under 
consideration is in the ballistic regime i.e. weakly disordered. Indeed, all the effects described 
here, result from the degeneracy, or quasi-degeneracy, of the ground state induced by the 
magnetic field in the case without Coulomb interaction. In principle, this always occurs 
for any cylindrical systems as carbon nanotubes, ring of carbon nanotubes or rings of usual 
semiconducting materials. 

The paper is organised as follows. In section II, the model without Coulomb interaction 
is presented and the origin of the crossing points induced by an applied magnetic field is 
described. In section III, we introduce the model with Coulomb interaction. In a first 
subsection, exact diagonalisation results are shown and the main interaction effects are 
discussed. In a second subsection, our variational ansatz is presented and some effects of 
the Coulomb interaction are shown to be well reproduced by our approximation. Last, in a 
third subsection, possible spin effects are analysed. 

II. NON INTERACTING ELECTRONS. ORBITAL EFFECTS 

The generic systems we consider along this work are rolled square lattices. In this 
section we start by neglecting the Coulomb interaction to focus on orbital effects only. The 
electrons are then described by the following nearest-neighbour tight-binding model where an 
uniform magnetic field, H, parallel to the cylindrical axis, is included via the Peierls-London 
substitution 



4 



N 2n M 

H = t ( C i+l, m ,<T C n, m ,*e 1 ^ + h.C.) + t (4 h rn+l,*C n ,m,a + h.C.). (1) 

n=l,(T=±5 m=l,(T=±5 

where is the magnetic flux through the section of the cylinder in units of the flux quantum 
0o (0o — hc/e). The two indices (ra, to), two integers, are the coordinates of the lattice sites: 
n is the coordinate along the circumference, 1 < n < N, and to the one along the cylindrical 
axis, 1 < to < M. The fermionic operator (c„ imiCr ) is the creation (destruction) 

operator of an electron at site (n, to) with spin a. The spectrum of the Hamiltonian (1), 
which depends continuously on the magnetic flux, is 

2lT 7T 

e M (0) = 2tcos(— (p + (/>)) + 2tcos(^— -g) (2) 

with p and g two integers such that —N/2 < p < N/2 — 1 and 1 < q < M. We have 
applied open boundary conditions at the ends of the cylinder. Note that the spectrum, and 
therefore every thermodynamic quantity, is periodic in flux, with periodicity 0=1 [3]. As 
the magnetic field is increased, the energy levels evolve and many level crossings appear 
[16] (Cf. Fig. 1). The details of the pattern of the crossing points is very complicated and 
depends on the geometry of the system and parameters of the model. 

At zero temperature, the ground state energy of the system with N e electrons is obtained 
by filling up successively the lowest energy levels according to the Pauli principle. Here, we 
restrict our study, almost exclusively, to the case of equal numbers of up and down spins i.e. 
S z = 0, where the lowest A^ e /2 levels are doubly occupied. The important variations of the 
energy levels with the magnetic field cause changes in the level occupation. As a consequence 
of that, at certain values of the magnetic field, C , a former excited state may become the 
new ground state. Such a switch of ground state are responsible for the appearance of cusps 
in the ground state energy curve. We can see an example in Fig. 2a for a cylinder with 
N = M = 10 and 80 electrons; the ground state energy as function of the magnetic flux 
shows 5 different cusps. At these particular points, the ground state is changed from a state 

> to a new one \ty > (i for 'in', and o for 'outf to do an analogy with scattering theory) 
which are both a Slater determinant built with the one-electron states, {(p P , q }, eigenf unctions 
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of the Hamiltonian (1). The two determinants differ only by the highest occupied level, fp iqi 
and <Pp qo (where the upperscript H is for Highest). 

One may see this orbital effect as a succession of scattering events where the time is 
replaced by the magnetic flux. The system of N e particles evolves freely until a particular 
'time', 0c, where the most energetic particle is scattered: <Pp iqi — > fp ,q - Depending on the 
one-electron state exchanged, one should distinguish between 

• 'forward scattering 1 (FS) where the two one-electron states, (pi,<&) and (p ,Q ), corre- 
spond to particles moving in the same direction along the circumference, 

• 'backward scattering' (BS) where the two one-electron states, (pi,qi) and (p ,Q ), cor- 
respond to particles moving in the opposite direction along the circumference. 

In general, a FS corresponds to smaller change in momentum, 5k = ([(pt — p )j[] 2 + [(<& ~ 
Qo) A^n"] 2 ) 1 ^ 2 ' than a BS. In Fig 2a, only the cusp at ~ 0.12 is associated with a FS event, 
all the others are BS events. 

The persistent current (PC) is a thermodynamic quantity given, at zero temperature, in 
terms of the ground state energy E by 

9E(4>) _ _ 2 ^ de pM ( 3 ) 



' PC 



The second equality arises for free electrons only and in the case where S z = 0. It can be 
shown that this derivative is proportional to the average of the current operator. This current 
yields an orbital magnetic moment which can be detected experimentally [4]. Obviously, 
the persistent current presents discontinuities for each value of the magnetic flux where the 
ground state energy shows cusps. Fig. 2b gives an example for the very same cylinder (N = 
M = 10, iV e = 80). Note that the persistent current in mesoscopic cylinders was studied 
with some details in the past. For instance, it was shown in [17] that its intensity strongly 
depends on the shape of the Fermi surface. This property has important consequences 
for carbon nanotubes [13]. Most of the other studies insist on the role played by disorder 
[15,18]. As a remark, one may say that the orbital effect described above, works to reduce 
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the persistent current: if one would keep the level occupation frozen and then, let evolve the 
ground state energy as function of the magnetic flux the resulting PC would be more than 
one order of magnitude larger. 

The discontinuities seen in the PC (Fig. 2b) is a general phenomenon, within the free 
electron picture, which affects any response function. For instance, in ref. [16], the static 
electric polarisability of diverse cylindrical systems was studied and shown to present these 
characteristics. It was then suggested to use this physical quantity to get some insights into 
the electronic structure of nanoscopic materials such as carbon nanotubes. 

In more realistic situations, the neighbouring energy levels (Fig. 1) are coupled by various 
interactions i.e. Coulomb interaction, disorder. In consequence, the levels will not cross each 
other but rather will come close and then repel in avoided crossings. Therefore, the different 
response functions of the system will not show discontinuities but rather abrupt changes at 
the position of the avoided crossings. Next we consider the case of long and short range 
Coulomb potential in the weakly interacting limit assuming the system without disorder. 

III. INTERACTING ELECTRONS. COULOMB AND SPIN EFFECTS. 

In this section, we consider the same cylinders pierced by a magnetic flux but, with 
interacting electrons. Therefore, we consider the following Hamiltonian 

H = H + H mt (4) 

where H describes free electrons in applied magnetic field and has been defined previously 
(see Eq. (1)). The second term introduces the Coulomb interaction 

Hint — — U( n ,m),(n' ,m')C n ,m,a C n' ,m' ,a' C n' ,m' ,a'Cn,m,a (5) 

(n,m)cr,(n' ,m')cr' 

where U(n,m),{n',m') is assumed to be either short range (Hubbard model) or long range. In 
the case of long range potential, we choose, to be specific, the Ohno potential known in the 
chemical literature where it is often used to describe n electrons in organic materials such 
as conjugated polymers or carbon nanotubes 
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U(n,m),{n>,m>) ~ U 2 ^ 

V 1 + a O r ( n ,m),(n',m') 

It includes a kind of effective screening due to the inner electrons via an effective screening 
length, a , with typical value of 0.61lA~ 2 [19]. T( n ,m),(n',m.') is the distance between two 
electrons siting in sites (n, m) and (n', m') given in angstrom; we choose the lattice units 
to be close to the usual carbon-carbon bond length in graphite, a = \AA. The functional 
form of the potential is then fixed and only U is kept as a variable. Apart the specificity 
of this particular potential, its main characteristics are rather general: a 1/r behaviour at 
large distances and an effective screening which prevents from any discontinuities at short 
distances. It could then also be taken as a reasonable Coulomb potential for other systems 
such as semiconducting rings or, for a very small (Hubbard model), metallic rings. 



A. Exact diagonalisation studies 

We have first considered very small cylinders for which it is possible to diagonalise exactly 
the Hamiltonian with interaction (Eq. (4)). Hereafter, particular results for a cylinder with 
N = 3, M = 2 and N e = 4 (two up and two down electrons) are discussed. For this filling, 
the ground state energy curve shows one crossing point in magnetic field in the case without 
interaction. The results are shown in Fig. 3 for increasing values of the Coulomb potential. 
Three effects of the Coulomb interaction can be seen. 

• First, as expected - but it is not clearly shown in the figure - the crossing point of 
the free-electron model is replaced by an avoided crossing: the Coulomb interaction 
opens a gap. Moreover, the avoided crossing becomes more and more pronounced with 
increasing values of the Coulomb interaction. 

• Second, the position in magnetic field of the (avoided) crossing point is shifted, C — > 
0c and the importance of the shift increases with the intensity of the interaction. In 
the following, we call this effect the 'Coulomb effect. 
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• Third, for large enough U, one notes the formation of a new plateau like structure in 
the ground state energy curve at the position of the (avoided) crossing point whose 
size increases with the interaction strength. Because of the appearance of this plateau, 
two new real crossing points are formed. We will see latter that this plateau may be 
explained by the spin degree of freedom: this is a 1 spin effect. 

These effects cause by the Coulomb interaction, will be studied in more details in the 
following two subsections using simple approximations allowing for the study of larger sys- 
tems. 

B. Two reference ansatz. Coulomb effect. 

We focus mainly on the weak-interacting regime where a mean field theory is supposed 
to be in principle, a good starting point. However, as we have already seen in the previous 
section, the free electron ground state becomes degenerate for some particular values of 
the magnetic field (cf. Fig. 2a). The Coulomb interaction contributes to mix the diverse 
electronic configurations thereby, at the vicinity of these points of degeneracy, one expects 
the ground state of the system to be a true many body state unable to be described by only 
one Slater determinant: a mean field description is then not appropriate. We propose, in 
this subsection, a simple variational ansatz going beyond a mean field treatment and able 
to capture some of the important many-body effects. 

It was suggested long ago [20] , to extent the usual Hartree-Fock theory for linear combi- 
nation of Slater determinants 

\^>=J2a k \^ k > (7) 
k 

where a k are variational parameters and \&k > are some chosen Slater determinants built 
from one-electron states determined by the variational principle [21,22]. The choice of Slater 
determinants entering the composition of \ip > is, of course, motivated by the problem under 
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studies. This kind of theory is particularly relevant in case of degeneracy such that appearing 
here for some values of the magnetic flux, C . 

Before studying the minimal ansatz of the form (7) relevant for our particular problem, 
we start by introducing new operators which are linear combinations of the site operators 
seen previously 

— a n,m C n,m,a (8) 
n,m 

where a l nm are complex coefficients. At best, these coefficients are determined after a vari- 
ational procedure as it will be done latter. In the free electron model, these operators 
correspond to the molecular orbitals, eigenfunctions of H . In any case, they are listed by 
increasing value of their corresponding energy. 

In the non-interacting case, the ground state energy shows points of degeneracy as func- 
tion of the magnetic flux. Let us first, analyse in details these different crossing points, this 
will help us to define our ansatz. For S z = 0, each of these points shows, at least, a fourfold 
degeneracy. The degenerate states are listed in the following. 

N 

i^>=n44i°> ( 9 ) 

i=i 

where N = N e /2 and |0 > is the vacuum without electron. This state is the state \ipi > 
mentioned above. The three other degenerate states are 

1^// >= A^A^A^A^ > (10) 

this is the state \ip Q >, 

\i>m>= A^ +n A m \^> (11) 

and 

Wiv >= A^ +n A m \^j > . (12) 

This analysis is illustrated in Fig. 4 where the four Slater determinants are represented 
schematically together with the behaviour of their corresponding energies as function of the 
magnetic field. 
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The electronic configurations > and \ipn > are closed-shell. With Coulomb interac- 
tion, they are expected to give the main components of the ground state at values of the 
magnetic flux sufficiently lower or higher than <f> c , respectively. The two other configura- 
tions, \ipni > and \ipiv >, are mono-excitations with respect to \ipi > or \ipu >. They are 
expected to play a role in the direct vicinity of C (Fig. 4). As a first trial, we consider 
only the closed-shell determinants in the linear combination (7); we will see, at the end, 
that this simple ansatz gives good results not too close to the degeneracies and is sufficient 
to describe (i) the avoided crossing formation and (ii) the Coulomb effect but not the spin 
effect. Within this approximation, we consider the following two- reference ansatz 

JV 

^ >= a \ijj > +P\i> n >= [al + pA^ +n A^A Nl A m ] J[ A^AQ > (13) 

i=i 

where 1 is the unit operator. The expansion coefficients, a and (3, and the set of coeffi- 
cients {a^ m } used to define the one-body wave functions (8) are now determined from the 
variational principle. More precisely, we look for an extremum of the ground state energy: 

'<«■>.{<->)= ^5^- <"> 

We use a two-step iterative procedure to determine the best wave-function \ip >. In 
a first step, the expansion coefficients a and f3 are determined by diagonalisation of the 
following two by two matrix 



1 <^\H\^r> <^j\H\^ n > 



(15) 



In a second step, the coefficients a l nm are determined keeping the expansion coefficients, 
a and (3 fixed. Minimisation of £ (a, /3, {a^ m }) with respect to the set {a l nm } gives the 
so-called generalised Brillouin theorem for multiconfigurational Hartree-Fock theory [21,22], 
valid for any multireference states (7) 

<^\[AtA s ,H}\^>=0. (16) 
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The Brillouin theorem for usual HF theory states that the matrix elements of the total 
Hamiltonian between the ground state and any mono-excitations cancel out. Its generali- 
sation to multireference HF theory is less clear: it means that the matrix elements of the 
total Hamiltonian between the ground state and some linear combinations of excited con- 
figurations cancel out. Using this condition, it is possible to calculate the set of coefficients 
{ a n, m } as we wn l see below. Then, the process is continued by repeating these two steps 
until convergence is reached. 

The first step is straightforward. For the second step, to proceed starting from Eq. (16), 
we follow closely the Ref. [21]. In order to make use of the generalised Brillouin theorem, it 
is convenient to rewrite the Hamiltonian as 

H = F — V + H int (17) 

where V is an effective one-body operator and, F is the Fock operator which we require to 
be diagonal 

F = H + V = ]T 8 rs e r Al ja A r , a . (18) 

r,s,a 

With this rewriting of the Hamiltonian, the generalised Brillouin theorem is used to deter- 
mine V and, thereby, to build the Fock operator. By using a multireference ansatz such as 
(7) and (13), the Fock operator takes a block structure 











F 

- 1 cc 


F 

A CO 


F 

- 1 ce 




F 

± CO 


F 

± oo 


F 

± oe 




F 

- 1 ce 


F 

- 1 oe 


F 

- 1 ee 



\ 



Let us name £ r Xo-, with 1 < r < N.M and a = ±1/2, the variational one-electron wave 
functions associated with the fermionic operators Al a and A T ^. £ r is the orbital part and 
Xa the spin part. They are listed by increasing values of energy. The one-electron operator 
F cc is defined in the subspace spanned by £ r with 1 < r < N — 1 (N — N e /2), F OQ in the 
subspace spanned by and £tv+i an d F ee in the subspace spanned by £ r with r > N + 1. 
It is important to stress that the effective Fock operator is non-uniquely defined by the 
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generalised Brillouin condition. Indeed, only the off-diagonal blocks are determined by 
using Eqs. (16) and (18). Additional assumptions are needed to fix the three diagonal 
blocks: here, we calculate these blocks as matrix elements of the one-particle operator F ce , 
procedure particularly appropriate for closed shell configurations [21]. 

In order to test the quality of the variational wave function (13), we have first considered 
very small cylinders and compare the approximate results with the results obtained by exact 
diagonalisation. Fig. 5 shows a comparison between the variational calculation and the exact 
result for weak interaction (U = O.lt and U = 0.2t). A good agreement is obtained which 
becomes excellent away from the crossing points. Indeed, Fig. 5 shows only the direct 
vicinity of the crossing point. As we can see, our simple variational ansatz (13) describes 
well the Coulomb effect i.e. the shift of the crossing point, but is less accurate at the 
crossing itself. The reasons for that could be the spin effect and/or the two open-shell Slater 
determinants not included in the ansatz (13). 

Before going further, let us discussed in some details the origin of the Coulomb effect 
well described by our ansatz. For simplicity, we neglect the off-diagonal term of the two by 
two matrix Eq. (15) - responsible for the opening of the gaps - and the Fock term keeping 
only the most important contribution of the Coulomb interaction, namely the Hartree term. 
As a last approximation, instead of considering the self-consistent procedure just described 
we keep only the first order correction. Without interaction, at every cusp, C , the system 
changes of ground state from |^ > to \$> >. The energies of these two states are Ei{4>) 
and E a ((f>). They are both two parabola-like curves that cross at <p c , Ei(<f) c ) = E Q ((f> c ). The 
electron densities of these two states are Pi(n, m, <f>) and p (n, m, 0), respectively, which differ 
in general i.e. Pi(n, m,(f>) ^ p (n,m, 0), 

p a (n,m,0) = 2 \ l P P c„q a ( n , m , ( t ) )\ 2 (20) 

OCC. 

where a = i/o. The first order Hartree corrections, Wi{4>) and W {4>), are given by 

W a {4>) = g H Pa(w,m,^)C/'( n>m ) ) ( n / >ro /)p Q (n / ,m',0) (l - ^5( n ,m),(n',m')) • ( 21 ) 

(n,m) ^ 
(n' ,m') 
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Since the densities differ, the magnitude of the Hartree corrections will also be different - 
in particular, at the crossing point, Wi(<f> c ) ^ W ((f) c )). In consequence of that, the two 
parabola, Ei(<f>) and E (4>), are shifted more or less strongly with Coulomb interaction to 
higher energies by the Hartree corrections, Wi(4>) and W {4>). Since the Hartree terms have 
different magnitude, the two-shifted parabola will cross at a different value of the magnetic 
flux, 0c, such that Ei(<f> c ) + Wi(<t> c ) = E ((j> c ) + W (<f> c ). This is the Coulomb effect illustrated 
schematically on Fig. 6 and observed, for instance, in the exact results shown in Fig 3 and 
Fig. 5. 

The magnitude of the shift is determined mainly by the difference between the Hartree 
contributions, 8W{4>) = Wi(<f>) — W (<j>), caused by the highest occupied states, f^ iqj an d 
<Pp q - If one uses more compact notation replacing the two pairs of indices, (pi,qi) and 
(p , Qo), by ki and k Q , respectively, and if one works in the momentum space instead of the 
real space, the mean contribution to 5W((f)) is given, at first order in perturbation theory, 
by the following equation 

5W{4>) oc £(KlXK+ q | 2 - |<| 2 t/ q |vW) ~ ECKl^K+ql 2 - Kr^lvW) 

q q 

(22) 

where the summation is done for states that are occupied in the ground state, U q is the 
Fourier transform of the Coulomb potential and Sq = k — ki. On the basis of this approxi- 
mate formula, one can draw some conclusions. First, one sees that the magnitude of 5W(4>) 
would be enhanced by a long range Coulomb potential which gives diverging contribution 
at small momenta, q, (for instance, the Coulomb interaction behaves as 1/q for a two di- 
mensional system) on the contrary to short range interaction that provides only constant 
contribution (independent of q). Second, If Sq is small, as it happens for Forward Scattering 
(FS), one may neglect in the formula (22) the difference between the one-particle states 
(p^. and ipf^. Then, SW(4>) becomes roughly proportional to the derivative of the Coulomb 
potential, dU q /dq, which is a diverging quantity at small values of q (for infinite lattice) 
in the case of long range potential. To conclude, (i) we expect a large Coulomb effect for 

14 



long range potential and, comparatively, no effect for short range potential, (ii) in the case 
of long range potential, since the difference between ki and k Q are smaller for FS crossing 
points, than BS crossing points, the Coulomb effect are expected to be stronger for FS than 
BS. These two conclusions are confirmed both by exact diagonalisation for small cylinders 
(see Figs 3, 11 and 12) and by using our variational ansatz as we will see below. 

We consider a bigger cylinder with N = 10, M = 10 and N e = 80 (40 electrons up and 
40 electrons down). Then, it is not possible anymore to do exact calculation and one has to 
rely on the approximate result. Important shifts of all the crossing points are obtained (cf. 
Fig. 7). As it is expected from the qualitative argument above, the most important shift 
is obtained for FS type of crossing points: in our example a shift of more than 5.1O~ 2 0o is 
obtained which could correspond to significant value of magnetic field. Last, note that with 
on-site interaction only, no differences compared to the non-interacting case can be detected 
at the scale of the figure. This confirms that the magnitude of the shifts is controlled by 
the range of the Coulomb potential, in agreement with our qualitative argument based on 
Eq. (22). The screening of the Coulomb potential, induced by the proximity of a metallic 
electrode, for instance, will tend to decrease the magnitude of the Coulomb effect. 

The corresponding persistent current is shown in Fig. 8 compared to the one obtained 
for the very same cylinder but without Coulomb interaction. The discontinuities appear 
shifted to higher magnetic flux due to the Coulomb effect. This is particularly sensible for 
the FS event for which, moreover, the magnitude of the discontinuity is enhanced. 

C. Spin effect 

It is clear from the results of exact diagonalisation (cf. Fig. 3), that the ansatz (13) 
is not sufficient for relatively high Coulomb interaction. According to the analysis of the 
crossing points, we have missed in the linear combination (7), at least, the two mono- 
excitations \ipin > and \ipiv >• In the non-interacting case (U = 0), these two configurations 
are degenerate in energy, but, in the interacting case (U ^ 0), they are mixed into two 
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components that split from each other: a singlet state, at high energy, and a triplet state, 
at low energy, with a gap of 2r, T being the exchange energy 

r= | < rl) in \H\il>iv > I- (23) 

The corresponding wave functions are given by the symmetric and antisymmetric linear 
combinations 

\i>T/s >= -j^M > T\i>iv >] = ^=[4r+i,T^,T T A n+i,i A n,iMi > ( 24 ) 

where the indices T is for Triplet and S for Singlet. The triplet state may be lower in energy 
than the variational ground state found in the previous subsection. 

In principle, a first improvement of the ansatz (13) would consist to add the components 
(11) and (12) to the trial wave function. This would increase the number of variational pa- 
rameters by two and, more importantly, since \ipni > and \ipiv > are both mono-excitations 
i.e. open shell configurations, this would force us to change the procedure of calculation. To 
give an analogy with usual HF theory, one would have to change from a restricted Hartree- 
Fock to unrestricted Hartree-Fock type of calculation [21]. We postpone this treatment to 
forthcoming work and, instead, adopt a kind of minimal strategy to get a first understanding 
of what could be the role of the spin at the vicinity of the crossing points. Starting from 
the variational solution (13), we build the triplet wave function (24) and then estimate its 
energy 

£ T (a,0,{a\ hm })=<il>r\H\il>r> (25) 

keeping the parameters a, ft and {a l nm } frozen. The results obtained for the triplet state 
are then not variational, but should give, nevertheless, useful indications about the total 
spin of the ground state as function of the applied magnetic field. At the vicinity of the 
crossing points, the triplet state can be lower in energy suggesting a possible spin transition 
induced by a longitudinal magnetic field. Such an example is shown in Fig. 9 for a very 
small cylinder with U = t. The two reference ansatz (13) reproduces well the shifting of 
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the crossing point (see also Fig. 5), but the additional 'plateau' induced by the interaction 
is missed with this approximation. However, it is partially described if one considers also 
the energy (25): on a certain interval of magnetic flux - corresponding to this plateau - the 
triplet state (24) is lower in energy. Our procedure is in a good agreement with the exact 
results and, therefore, strongly suggests the following sequence for the total electronic spin 
at the vicinity of each point of (near) degeneracy and for increasing values of magnetic field, 
S — — — 0. At each spin transition appears a new real crossing points. 

Fig. 10 shows the corresponding persistent current compared with the one obtained for the 
non-interacting case. The long-range Coulomb interaction shifts the discontinuities - this is 
the Coulomb effect - and is responsible for the appearance of a new 'plateau' - the Triplet 
plateau - which may be detected experimentally. 

The above result may be explained by an usual qualitative argument based on the Hund's 
rule. In case of double degeneracy, or near degeneracy, such as encountered at the direct 
vicinity of each crossing point, and for S z = 0, the Coulomb interaction will favour the 
state where the two degenerate levels, <Pp iqi an d ¥p , qo -> are mono-occupied. Indeed, it costs 
more Coulomb energy to doubly occupied one of these two levels. Moreover, the repulsive 
interaction will favour the state with maximum total spin, here S = 1, since the coordinate 
wave function is then anti-symmetrised which allows to gain the exchange energy. The two 
highest electrons will occupy the levels <Pp uqi and <fp qo with total spin one, as long as the 
energy difference between the levels remains smaller than the energy gain due to Coulomb 
interaction. At this point, it is important to remember that we have not considered the 
Zeeman interaction which may change the picture at high magnetic field. 

Before doing the same analysis for a bigger cylinder, we present some exact diagonal- 
isation studies of the very same small cylinder but with a short-range Coulomb potential 
i.e. namely by using the Hubbard model. The results are summarised in Fig. 11 where 
the ground state energy is shown as function of the magnetic flux for several values of U . 
These results have to be compared with the ones of Fig. 3 presenting the same quantity but 
using the long-range Coulomb potential (Eq. 6). One sees, first, that the Coulomb effect 
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is considerably reduced by the use of an on-site term only: there is almost no shifting of 
the crossing points. Second, the spin effect is, on the contrary, sensibly enhanced in the 
sense that (i) the Triplet plateau appears for lower values of U than in the case of long 
range Coulomb potential and, (ii) this plateau can be much more extended in magnetic 
flux. It seems that there is a kind of competition between the two effects. The more the 
Coulomb effect is important the more the size of the Triplet plateau is reduced. This can 
be understood by the following qualitative argument. 

The difference in energy between the variational ground state and the Triplet state is 
given by 

F(a, [3, {< m }) = £(a, (3, {< m }) - &r(a, (3, {< m }) = 

(26) 

\ E(n,m),(n>, m >){pi{n, m) - p (n, 77l)) £/(„, m ),(n',m') (Pi « 77l') - p (n' , m')). 

If T > 0, the Triplet state is the ground state of the system. The on-site part of the Coulomb 
potential gives always a positive contribution to T in favour of a Triplet ground state. On 
the contrary, since the difference in density, 5p(n,m) = pi(n,m) — p (n,m), is in general 
an oscillatory function, the long range part of the Coulomb potential could give a negative 
contribution to T . Therefore, as it is seen in exact diagonalisation studies, the spin effect is 
expected to be more pronounced for screened interaction. 

Next, we have considered a bigger cylinder with N — 10, M — 10, N e = 80 and S z = 0, 
and evaluated the variational energy (Eq. (14)) and the corresponding triplet energy (Eq. 
(25)). For the long range Ohno potential (Eq. (6)), there is no range of magnetic flux where 
the Triplet state would be significantly favoured within our approximation - at least, for 
relatively low values of U (up to U — 2\t\ approximately). On the contrary, for the Hubbard 
potential, there is a Singlet — > Triplet transition at every crossing points for strong enough 
U. This is shown in Fig. 12 for U — \t\. One can notice, a large Triplet plateau for values 
of the magnetic flux corresponding to the FS point. 

Last, we have considered exclusively the case with equal numbers of up and down elec- 
trons up to now i.e. S z = 0. The conclusions will be different with an unpaired electron i.e. 
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S z = ±1/2. Indeed, in this case, at each crossing point the ground state of the free electron 
model is only twofold degenerate and not fourfold as it is the case for S z = 0. The spin 
is then not expected to play any role. Therefore, with Coulomb interaction, the avoided 
crossing and Coulomb effect phenomena should persist but, on the contrary, the spin effect 
described above should not. We have down exact calculations for very small cylinders. Fig. 
13 shows examples for N = 3, M = 2, U = 100|£|, 2 spin up and 1 spin down compared to 
the case with 4 electrons and S z = 0. One clearly sees that, for both long and short range 
Coulomb potential, the spin effect doesn't occur in the case with one unpaired electron. 

As a remark, one may say that the spin transition found at the vicinity of the crossing 
points reminds the fractional Persistent Current analysed for purely one-dimensional rings 
[14] . Without interaction, the ground state energy of these systems is periodic with periodic- 
ity 0o [16]. With infinite interaction, the spin-charge separation phenomenon occurs and the 
spin excitations of the system contribute to change drastically the behaviour of the ground 
state energy by changing its periodicity to <f>o/N e where N e is the number of electrons. It 
would be, of course, very interesting to be able to extrapolate what could happen for the 
multichannel systems studied here but in presence of very large interaction. 

Before the conclusion, one may add that similar spin sequences, as the one pointed out in 
this subsection, were already predicted to occur in short armchair carbon nanotubes under 
the influence of an inhomogeneous gate potential [23]. The carbon nanotubes have two 
families of levels which react differently to this applied potential, creating level crossings. 
Then, for even number of electrons, in the very same way as the situation described here, 
the first Hund's rule drives the system from a singlet (S = 0) to a triplet state (S = 1) 
at the vicinity of the crossing points. The same phenomenon was also described to happen 
in quantum dots where the signature of the spin transition was expected to be seen in the 
variation of the Coulomb-blockade peak positions [24]. 
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IV. DISCUSSION AND CONCLUSION 



We have considered in this work electrons moving on the surface of nanoscopic cylinders 
described as rolled square lattices. Without disorder - or with a weak disorder - a static 
magnetic field applied along the cylindrical axis induces points of degeneracy - or quasi- 
degeneracy. Therefore, with Coulomb interaction, the ground state of the system becomes, 
at the vicinity of these points, fundamentally a many body state unable to be described by 
any mean field theory based on only one Slater determinant (as usual Hartree-Fock theory). 
To circumvent this difficulty, we have studied a simple variational ansatz which is a linear 
combination of the two closed shell Slater determinants that become degenerate. On top of 
that, we have also considered the lowest triplet state built out of the variational solution. 
As a result, we have found three effects induced by the electron-electron interaction at the 
vicinity of these points of degeneracy. 

• Avoided crossings. The Coulomb potential induces couplings between the degenerate 
configurations. As a result, a gap is opened. 

• Coulomb effect. With repulsive interaction the diverse energy curves, corresponding 
to diverse Slater determinants, as function of the magnetic flux, are shifted up by 
the Hartree contributions. Because the magnitude of these contributions differs from 
configuration to configuration, the Coulomb interaction shifts the position in magnetic 
field of the (avoided) crossings. 

• Spin effect. When degeneracy occurs and for S z = 0, the Hund's rule can drive the 
system to a triplet state. Therefore, at the vicinity of the degeneracy points, the 
Coulomb interaction can induce a Singlet — > Triplet transition follows by a Triplet — > 
Singlet one. In the case of one unpaired electron, S z = ±1/2, the spin effect doesn't 
occur. 

The Coulomb effect is due mainly to the long range part of the Coulomb potential. On the 
contrary, the spin effect is favoured by short range interaction. Therefore, a screening of the 
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repulsive interaction induced, for instance, by a gate electrode could change the electronic 
properties of the system. 

The electron-electron interaction, on the one hand, mixes the different configurations 
and replaces the crossing points of a free electron theory by avoided crossings but, on the 
other hand, the total spin of the system may be changed due to the Hund's rule creating new 
real crossing points. Thereby, any response function, providing that the electronic structure 
of the cylinder is preserved, are expected to show an abrupt change at the position of these 
crossing points where the total spin is changed. This should be the case in measures of 
persistent current [4], but also, in the static electric magnetopolarisability studied in [16] 
and already measured for an ensemble of metallic rings [25] and, as a last example, in 
magnetoconductance measurements, with bad contacts to the electrodes, such as the one 
done in [26] where multiwall carbon nanotubes behave as quantum dots. Similar effects were 
studied in [24] for quantum dots, where the spin transitions were shown to give characteristic 
signatures in the Coulomb-blockade peak positions. 
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FIG. 1. Lowest energy levels of a cylinder with N = 100 (number of sites along the circumfer- 
ence) and M = 100 (number of sites along the cylindrical axis) as function of the magnetic flux 4>. 
Numerous level crossings appear at values of (j> which depend on the cylinder geometry. 
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FIG. 2. Cylinder with N = M = 10, N e = 80 and S z = 0. (a) Energy of the ground state 
as function of the magnetic flux. Doing an analogy with scattering events, the different cusps are 
classified into Forward (FS) and Backward (BS) Scattering types (see text). Here, only the cusp 
at cf) ~ 0.12 is a FS, the other ones are all BSs. (b) The corresponding persistent current. The 
discontinuities are more pronounced for BS than FS. 
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FIG. 3. Exact ground state energy as function of the magnetic flux of a small cylinder (N = 3, 
M = 2, N e = A and S z = 0) for increasing values of the long-range Coulomb potential, U, in 
units of \t\. The crossing point is shifted to higher magnetic flux with increasing interaction. For 
strong enough interaction, there is appearance of a new plateau that may be explained by spin 
excitations. The curves are shifted up for clarity 
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FIG. 4. Picture of a crossing point in the free electron case. For S x = 0, each crossing point is 
fourfold degenerate. The four Slater determinants, \ipi >, \ipu >, |V>m > and \ipiv > ( see text) 
are shown schematically together with the behaviour of their corresponding energies as function 
of the magnetic flux. The full rectangles represent all the doubly occupied levels except the two 
highest ones, ^Pp uqi and l Pp 0:qo , which are explicitely shown. In the representation of the Slater 
determinants, the energy levels are kept fixed at a particular value of </>, for simplicity. 
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FIG. 5. Zoom of the ground state energy of a small cylinder (N = 3, M = 2, N e = 4 and 
S z = 0) with long range Coulomb interaction in the direct vicinity of the crossing point. The 
full curve and the dashed curve are the exact results for U = 0.2 and U = 0.1 (in units of \t\), 
respectively. The full dots and the full squares are the results obtained with the two-reference 
variational ansatz (Eq. (13)) for U = 0.2 and U = 0.1 respectively: the Coulomb effect is well 
reproduced but not the appearance of a plateau which is the result of spin excitations not contained 
in the ansatz. The curves corresponding to U = 0.1 are shifted down for clarity. 
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FIG. 6. Schematic representation of the Coulomb effect. The dashed curves represent the 
ground state energy of the system without Coulomb interaction: there is a cusp at <f> c correspond- 
ing to a change of ground state \ipi >^ \ip >; at this point the ground state is degenerate, 
Ei{4> c ) = E (4> c ). The full curves represent the ground state energy with Coulomb interaction: the 
two initial parabola are shifted up by the Hartree contributions, Wi and W Q , which have different 
amplitudes; consequently, the two interacting ground states are exchanged at a different value of 
the magnetic flux, cf) c , such that Ei((f> c ) + Wi((f> c ) = E ((f> c ) + W (<fi c ). In other words, the long 
range Coulomb interaction produces a shift of the cusp, <p c — ► <f> c . 
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FIG. 7. Ground state energy as function of the magnetic flux for a cylinder with N = M = 10, 
iV e = 80 and S z = 0. The full curve is obtained with the two reference ansatz (Eq. (13)) for 
U = 0.08 in units of \t\, the dashed curve is for the case without interaction (U = 0). All the cusps 
are shifted to higher magnetic flux by the long-range Coulomb interaction; this is the Coulomb 
effect, particularly important in the FS case. One may notice that the same calculation, but 
with an on-site interaction only (Hubbard model), gives approximately the same result than the 
non-interacting one. 



30 



10 




-10 1 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 

FIG. 8. Persistent current as function of the magnetic flux for a cylinder with N = M = W, 
N e = 80 and S z = 0. The full curve is obtained with the two reference ansatz (Eq. (13)) for 
U = 0.08 in units of \t\, the dashed curve is for the case without interaction (U = 0). The 
discontinuities are all shifted to higher magnetic flux by the long-range Coulomb interaction. For 
the FS case, the shift is very pronounced and the importance of the discontinuity is enhanced 
compared to the non-interacting case. 
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FIG. 9. Zoom of the ground state energy of a small cylinder with TV = 3, M = 2, N e = 4 
and S z = 0, with long range Coulomb interaction, around the crossing point region. The full 
curve is the exact result for U = \t\. The open circles denote the result obtained with the two 
reference ansatz (Eq. (13)). The full dots denote the result obtained by taking the minimum 
between the values given by the two reference ansatz and the ones given by the triplet built out of 
the variational solution (Eq. (25)). This approximate results suggest transitions of the total spin, 
£ = 0— >5 = 1— >5 = 0. The Coulomb interaction, by changing the total spin, induces new real 
crossing points. 



32 



LU 
"O 



-10 




4> 



FIG. 10. Persistent current as function of the magnetic flux for a cylinder with N = 3, M = 2, 
iV e = 4 and S z = 0. The full curve is the exact result for U = \t\ and the dotted curve is for 
the non-interacting case. The long range Coulomb interaction induces a shift of the discontinuity 
plus appearance of a new plateau. This plateau can be explained by singlet-triplet transitions (see 
text). 
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FIG. 11. Ground state energy of a small cylinder with N = 3, M = 2, N e = 4 and S z = 0, with 
short range interaction (Hubbard model) for increasing values of U, in units of \t\, obtained with 
exact diagonalisation. Compared to the results with a long-range potential (Fig. 5), the Coulomb 
effect is strongly reduced but the spin effect is enhanced, in the sense that it becomes relevant for 
lower values of U and, that the Triplet plateau is larger. The curves are shifted down for clarity 
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FIG. 12. Ground state energy as function of the magnetic flux for a cylinder with N = M = 10, 
iV e = 80 and S z = with a short range electron-electron interaction (Hubbard model) and 
U = 0.4\t\. The full curve is the result given by the two reference ansatz (Eq. (13)), the dashed 
curve is the energy of the triplet state (Eq. (25)) built out of the variational ansatz. Around each 
avoided crossing, there is a sequence Singlet — ► Triplet — ► Singlet: for some intervals of magnetic 
flux, the system is in a triplet ground state. This interval is very large in the vicinity of the FS 
point. These changes of total spin give new real crossing points. On the one hand, the Coulomb 
interaction replaces the crossings of the free electron description by avoided crossings, but, on the 
other hand, creates new crossing points due to spin transition. A zoom on the last crossing points 
is shown in the inset. 
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FIG. 13. Ground state energy of a small cylinder with N = 3, M = 2 and U = 100|i|. The 
dotted and dashed curves are for iV e = 4, S z = 0, short and long range Coulomb potential, 
respectively. The dot-dashed and full curves are for N e = 3, S z = 1/2, short and long range 
Coulomb potential, respectively. The spin effect doesn't occur in the case with one unpaired 
electron (see text). 
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